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Bm 

EVr 
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F. 

Fj 
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5,7, 

Vn 
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L 

L. 

Lx,Lj 

Nr 

fc*> 

P 
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Constants 

Membrane  number,  see  Eq.  {16) 

Bond  number,  see  Eq.  {16) 

Amplitude  of  the  n*  sloshing  mode 
Potential  of  the  solid-li quid-gas  interface 
Potential  of  the  edge  load  acting  on  the  rim  of  the 
tank  bottom 

Gravitational  acceleration,  time  dependent,  see  Eq.  (6a) 
Mean  local  gravitational  acceleration 
Amplitude  of  the  imposed  axial  acceleration 
Nondimenaional  gravitational  acceleration,  see  Eq.  (15) 
Membrane  thickness 

Nondimensional  free  surface  shape,  see  Eq.  (15) 
Functionals,  see  Eqs.  (36) 

Bessel  functions  of  first  kind 

th  - 

n  root  of  the  equation 

Depth  of  liquid 

Nondimensional  depth  of  liquid,  see  Eq.  (15) 

Pressure  energy  in  nondimensional  form,  see  Eq.  (41) 
Lagrangians  in  nondimensional  form,  see  Eqs.  (43),  (44) 
Matrix 

Midplane  stress  resultant 

nt!l  sloshing  frequency  for  rigid  tank,  see  Eq.  (57) 
Pressure 

Nondimensional  pressure,  see  Eq.  (15) 
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Radius  of  the  tank 

Cylindrical  coordinates,  see  Fig.  1 
Nondimensional  cylindrical  coordinates,  see  Fig.  1 
Surfaces,  see  Fig.  1 
Time 

Modal  matrix 

Transverse  deflection  of  membrane 

Nondimensional  transverse  deflection  of  membrane,  see  Eq.  (15l 

Vertical  coordinate,  see  Fig.  1 

Nondimensional  vertical  coordinate,  see  Eq.(15) 

Nondimensional  amplitude  of  the  imposed  axial  acceleration, 
see  Eq.  (15} 

Constants,  see  Eq.  (59) 

Boundary  curves  of  Sj,  S3  respectively 
Constants,  see  Eq.  (66) 

Free  surface  shape 
Asimuthal  coordinate 
Mass  ratio,  see  Eq.  (16) 

Constants,  see  Eq.  (54a) 

Constants,  see  Eqs.  (51),  (52) 

Density  of  membrane 
Density  of  liquid 
Surface  tension 

Nondimensional  time,  see  Eq.  (15) 

Velocity  potential 

Nondimensional  velocity  potential,  aee  Eq.  (15) 

Forcing  frequency 


ABSTRACT 


The  stability  of  a  fluid  contained  in  a  circular  cylindrical 
tank  with  a  flat,  flexible  bottom  under  a  periodic  axial  excitation 
ia  studied.  An  analytical  difficulty  for  the  solution  of  the  linearized 
equations  in  the  form  of  infinite  aeries  is  discussed.  A  variational 
approach  is  formulated.  An  approximate  solution  results  in  a  pair 
of  coupled  ordinary  differential  equations  with  periodic  coefficients. 
A  method  of  handling  the  stability  of  the  solutions  of  such  a  system 
of  equation  is  presented.  Numerical  results  will  be  discussed  in  a 
later  article. 
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Eigenvalues,  see  Eq.  (5?) 

Frequency  parameter  for  the  membrane,  see  Eq.  (16) 
Frequency  parameter  for  surface  tension,  see  Eq.  (16) 
Two  dimensional  Laplace  operator  = 

in  cylindrical  coordinates  ^ 

Three  dimensional  gradient  operator 

in  cylindrical  coordinates.  For  the  free  surface  and 

the  membrane,  ignore  the  term.} 
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1.  Introduction 

Dynamics  of  large  liquid-fuel  rockets  naturally  involves  the 
motion  of  a  liquid  in  a  flexible  container.  The  symmetric  modes 
of  the  fluid  motion,  which  influences  the  pressure  at  the  tank  bottom,  -  , 
and  thence  influencing  the  pressure  in  the  pump,  combustion  chamber, 
and  thrust  and  rocket  acceleration,  have  an  Important  effect  on  the 
rocket-structural  dynamics.  In  some  instances  the  longitudinal 
oscillations  were  so  serious  as  to  affect  the  safety  of  the  vehicle. 

For  this  reason  the  analysis  of  the  forced  oscillations  of  the  liquid 
container  is  important. 

At  ground  level  perhaps  the  effects  of  the  flexibility  of  the 
tank  wall  and  the  surface  tension  of  the  free  surface  are  negligible 
on  fuel  sloshing.  At  rednced  gravity  conditions  these  effects  will 
become  more  evident.  It  is  the  purpose  of  this  article  to  evaluate 
the  effects  of  tank  flexibility  and  surface  tension  on  the  stability  of 
liquid  motion  in  the  symmetric  modes. 

Sloshing  of  liquids  has  been  studied  by  many  authors.  Most 
of  them  considered  rigid  containers,  see  Ref.  4.  However,  Miles 
(Ref.  1)  considered  bending  modes  of  a  flexible  container.  Blelch 
(Ref.  2)  investigated  the  longitudinal  modes  approximately.  Recently  ^ 

Bhuta  and  Koval  (Ref.  3  and  15)  studied  the  coupled  oscillations  of  a  t 
liquid  in  a  tank  with  a  flexible  bottom.  They  defined  the  normal  modes 
of  the  system,  and  treated  the  orthogonality  and  expansion  theorems. 
However,  the  difficulty  concerning  the  convergence  of  several  infinite 
series,  to  be  explained  in  Section  5  below,  was  not  considered.  In 
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Refs.  16  and  17,  Bhuta  and  Teh  considered  the  problem  of  arbitrarily 
assigned  velocity  distribution  on  the  tank  bottom;  in  this  case  the 
difficulty  referred  to  above  does  not  appear. 

On  the  other  hand,  there  is  substantial  literature  about  the 
Influence  if  surface  tension  on  sloshing;  see  Teh's  bibliography.  Ref.  18, 
and  papers  by  Bond  and  Newton  (Ref.  19),  and  Reynolds  (Ref.  20).  How¬ 
ever,  most  of  these  studies  are  concerned  with  free  oscillations,  very 
little  has  been  done  about  the  Influence  of  surface  tension  on  forced  j 
oscillations.  No  work  seems  to  have  been  done  on  the  coupling  with 
the  flexibility  of  the  tank. 

In  the  present  paper  a  circular  tank  with  a  flexible  bottom  under 

w* 

vertical  periodial  excitation  is  studied.  The  problem  was  first  formulated 

i 

in  the  form  of  differential  equations  and  then  in  the  form  of  a  variational 
principle.  An  approximate  solntion  is  presented,  which  results  in  a 
pair  of  coupled  ordinary  differential  equations  with  periodic  coefficients. 
The  stability  of  the  solutions  of  these  equations  is  discussed. 
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2.  Statement  of  the  Problem 

A  circular  cylindrical  container  with  rigid  side  wall3  and  a 
Oat.  flexible  bottom  contains  a  liquid  with  a  free  surface.  The  tank 
walls  are  subjected  to  an  oscillatory  axial  acceleration,  in  addition 
to  a  constant  mean -local-gravitational  acceleration  which  is  directed 
along  the  axis  of  the  cylinder.  A  gas  with  constant  pressure  exists 
i«  .1  above  the  liquid  surface.  No  external  force  acts  underneath  the  tank 
bottom.  The  situation  is  pictured  in  Fig.  1.  The  problem  is  to 
determine  the  motion  of  the  liquid;  in  particular,  its  stability. 

The  fluid  properties  including  the  surface  tension  are  assumed 
to  be  uniform,  constant,  incompressible,  and  imriscid. 

The  mean  free  surface  of  the  liquid  shall  be  assumed  to  be  a 
plane  perpendicular  to  the  cylinder  axis.  In  low  gravity  and  finite 
surface  tension  one  may  have  to  consider  a  curved  mean  free  surface. 
The  governing  criterion  is  the  Bond  number  defined  below.  In  thiB 
paper  we  shall  assume  that  the  Bond  number  is  sufficiently  largo  so 
that  the  free  surface  is  approximately  a  plane.  The  case  of  low  Bond 
number  will  be  investigated  in  a  separate  article. 

As  a  further  simplification  we  assume  that  the  deviation  from 
the  static  equilibrium  condition  is  small,  so  that  the  deflections  of 
the  free  surface  and  of  the  tank  bottom,  the  fluid  velocity,  and  hence 
the  velocity  potential,  may  be  considered  as  infinitesimal  quantities 
of  the  first  order.  Under  this  assumption  all  the  equations  can  be 
linearized,  and  the  mathematical  problem  is  relatively  simple.  A 
number  of  Interesting  nonlinear  problems  are  ruled  out  by  this  assump¬ 
tion.  But  as  an  investigation  of  the  initial  tendency  toward  instability, 
the  linearized  theory  should  be  adequate. 
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3.  Mathematical  F orm’dlation  •' 

- - * - 

Consider  a  quantity  of  inviscid  liquid  situated  in  a  cylindrical 
container  of  radius  Tq  as  is  shown  in  Fig,  1.  The  cylindrical  polar 
coordinate  system  is  chosen  so  that  the  positive  z -direction  is  directed 
upward  away  from  the  liquid,  the  zero  on  this  axis  being  fixed  on  the 
mean  free  surface.  On  assuming  the  fluid  to  be  inviscid,  incompres¬ 
sible,  and  the  motion  to  be  irrotational,  the  equation  of  continuity 
may  be  expressed  in  terms  of  the  velocity  potential  <}>  : 


and  the  velocity  components  u,  \!  ,  w  are 

The  usual  subscript  notation  is  used  to  denote  partial  differentiation. 

The  kinematic  conditions  at  the  tank  walls  and  the  free  sur¬ 
face  are 


-  -a<b 

on 

r^r0 

(3) 

vs5i-^ 

ag-ar 

&n 

(4) 

vu-fU  4 

cr»v 

(5) 

whore  w  denotes  the  deflection  of  the  tank  bottom,  and  denotes  the 
deflection  of  the  free  surface,  both  positive  in  the  positive  z-direction, 
and  both  are  assumed  to  be  infinitesimal. 

Since  the  motion  is  irrotational,  Bernoulli's  equation  is  satis¬ 
fied  throughout  the  liquid  domain.  In  particular,  at  the  free  surface. 


we  have 
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\  -  -  4  (y40"-  ^  r?  -  ^  ctt) 


ft  a  -  i  •  -  —  '  (6) 

gltl- -+<3,  OO  cot  (6a) 

where  c(t)  is  an  arbitrary  function  of  time,  gg  is  the  mean  local  gravi¬ 
tational  acceleration,  g  j  coscot  is  the  imposed  axial  acceleration?  both 
taken  as  positive  if  they  are  directed  toward  the  tank  bottom  (along  the 
negative  z-directicn),  and  p  is  the  pressure  just  inside  of  the  interface. 
The  pressure  p  is  related  to  the  pressure  just  outside  of  the  liquid,  p^., 
by  the  relation 


4  -  \>  =  <r  X  (7) 

where  (j"  is  the  surface  tension,  and  His  the  total  curvature  of  the 
free  surface.  In  linearized  form,  under  the  assumptions  that  VrQ  <<.  1 
and  jgrad  l|  |  1,  we  have 


If  the  pressure  of  the  gas  p^,  is  a  constant,  then  without  loss  of  gener¬ 
ality  we  may  set  p^.  =  0.  The  function  c(t)  can  be  absorbed  in  We 

can  also  neglect  l^fin  Eq.  (6)  and  evaluate  4*  j  on  the  surface  z  =  0 
under  the  scheme  of  linearization.  Thus  we  obtain  the  linearized  free 
surface  condition, 


-Sir)1.,  — 


r  ? - -  a*  7 


Here  we  just  write  out  a  special  form  of  imposed  axial  acceleration. 
The  method  developed  below  can  be  applied  to  a  general  periodic  im¬ 
posed  axial  acceleration. 
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Simiiarly,  Bernoulli's  equation  gives  the  pressure  on  top  of 
the  bottom  wall 


No  other  forces  will  be  assume  d  to  be  acting  on  the  tank  bottom.  If 
the  tank  bottom  is  very  thin  and  is  prestressed  so  that  it  behaves 
like  a  membrane,  then  the  equation  of  motion  of  the  bottom  is 

urvAw=  +  ,  (in 

where  Nf  is  the  tensile  stress  resultant  in  the  tank  bottom,  and  is 
assumed  to  be  a  constant,  f  is  the  density  of  the  tank  bottom 
material,  h  is  the  tank  bottom  wall  thickness,  so  that  |3h  is  the 
mass  per  unit  area  of  the  tank  bottom. 

A  combination  oi  (10)  and  (11)  gives  the  linearized  equation 
of  motion  of  the  elastic  Dottom  as  a  membrane 


In  reality,  a  tank  with  flat  bottom  will  develop  both  bending 
and  stretching  stresses  under  fluid  pressure.  Eq,  (12)  is  a  good 
approximation  only  if  a  membrane  tension  is  built-in  at  the  edges  by 
stretching  the  bottom  onto  a  rigid  cylinder  before  the  two  are 
welded  together. 
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It  is  necessary  to  specify  the  boundary  conditions  for  and  w 
at  the  edge  r  =  .  We  choose 

O  whea  r  ~  r*  (13> 

~b  _  0  when  ^  _  y-  (14) 

The  last  condition  is  a  special  case  of  zero  capillary-hysteresis. 

It  is  consistent  with  the  simplifying  assumption  that  the  undisturbed 
free  surface  is  a  plane  z  *  0.  In  very  low  gravity  condition  the  mean 
free  surface  is  curved  and  Eq.  (14)  should  be  replaced  by  the  condition 
9^/Sr  *  at  the  wall,  where  y  is  a  physical  constant. 

These  equations  define  the  linear,  inviscid  problem  of  sloshing 
under  appropriate  initial  or  periodicity  conditions. 
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4.  Dimensionless  Equations 

Taking  the  radius  of  the  cylinder  T»  as  the  characteristic 
length,  the  gravitational  acceleration  as  the  characteristic  accel¬ 
eration,  and  oo  as  th'  '-haracteristic  frequency,  we  define  the 
dimensionless  variables 


R=  JL 


rv  ,  2=^  , 

5  -  sr  H=  ^  ,  V\l- 


»  -  5. 

71  , 


w 

To 


^fCc) 


_  a.-r  a.taot 


a 


—  r>- 

•  >  '  f :o>v  > 


(15) 


the  dimensionless  parameters 
Bond  number  =  •=.  — ■ 

P* 

Membrane  number  a  P,  _  =  — *-■  *- 
M  Nr 

X.  p  *J" 

Frequency  parameter  for  surface  tension  =  _Q_^.  —  —-*■■■* - 

Frequency  parameter  for  the  membrane  =  —  ft  V?- 

Nr 


(16) 


Mass  ratio  =  A.  — 
and  the  operator 

Then  the  equations  become 

=  0 

V"H  -  EVaH  =  o 


(17) 


(18) 


(19) 
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jE^toVvj  = O 


(20) 


with  the  boundary  conditions 

i£-_  Q 
zR  ° 

2>W 

T2=  ax: 

ai  aH 
a?"  at 

W  —  o 


a 

a 


-  o  °a 


on 

(21) 

on 

?=-L 

(22) 

on 

(23) 

on 

R-  \ 

(24) 

on 

\  (assuming  ^=>  0) 

(25) 

Eqa.  (18)  -  (25)  show  that  the  problem  of  sloshing  depends  on 
the  parameters 

j  2><r,  crt  )  >  X  ,  l_ 

Those  dimensionless  parameters  are  not  all  independent;  for 


Hence 


=  2k 

2* 


(26) 


(27) 


However,  we  retain  the  symbols  j  'H.H,&KlbecaU8e  these  two 

pairs  of  parameters  are  not  likely  to  be  both  important. 
t  the  tank  bottom  is  rigid 

»f  the  surface  tension  has  no  effect. 
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5.  rtn  Analytical  Difficulty 

It  is  not  easy  to  construct  an  exact  solution  for  the  system  of 
equation!  (18)  -  (25).  To  construct  a  solution  in  the  form  of  an  infinite 
senes,  as  a  direct  extension  of  Benjamin  and  Ursell's  solution  for 
a  rigid  tank,  encounters  certain  basic  difficulties.  Consider  symmetric 
modes  of  motion  in  which  are  independent  of  the  angular 

coordinates  6  .  A  solution  of  (18)  may  be  posed  as 

1=  v  w 

Then  Eqs.  (22)  ana  (23)  give 

|2„ 

W  =  I  aco]  JsM)  <3°) 

Both  Eqs.  (21)  and  (25)  are  satisfied  if  the  are  the  roots 

of  the  equation 


n=  I,  2(  i . 

£q.  (24)  is  satisfied  by  taking 


(31) 


djot  J.  [dj(r)  -  ChCC^  J.IM  =  0 


(32) 


ffR)  is  the  static  deflection  of  the  membrane.  We  assume  that  if  B. , 

M 

i«  positive,  Xbm  is  less  than  the  first  root  of  JQ(x)  »  0,  namely, 

2.4048.  To  satisfy  Eqs.  (19),  (201,  we  substitute  $  ,  H,  W  from 
(28)  -  (30),  collect  terms  and  represent  the  left  hand  side  as 
Fourier-Bessel  series  in  J  (  ‘t  R).  Since  the  series  vanish. 
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every  coefficient  if  XO^R),,  _ ,  must  vanish.  Thus  from 

(19)  we  obtain  the  necessary  conditions 

Jit  CJttH-  E>,  W  <3st«  =  Q  (33a) 

jit  0.  (33b) 

And  from  (20)  and  (33a) 

Sll  ice)  (34a) 

+ IMbl 

Thus  all  differential  equations  and  boundary  conditions  are  satisfied 
by  the  assumed  form  of  |-^  provided  that  C,  ,  d„  satisfy 

Eqs.  (32),  (33),  and  (34). 

A  difficulty  become  s  apparent  when  one  examinos  the  Eqs.  (32)- 
(34).  If  we  truncate  the  infinite  series  (28)  -  (30)  to  Kj  terms,  we  see 
that  Eqs.  (32)  -  (34)  always  impose  N  +  l  conditions  on  k)  unknowns, 
which  in  general  have  no  solution. 

The  difficulty  appears  to  have  risen  from  the  condition  VJta  O 
at  the  edge  of  the  tank  bottom.  The  same  difficulty  would  have 
appeared  in  Benjamin  and  Ur  sell1  s  rigid  tank  case  if  one  assumes  a 
"stuck"  condition  at  the  edge  of  the  free  suriaco,U=o  aksitR**, (physically 
realizable  with  certain  fluid  and  wall  material).  Mathematically,  the 
condition  W=o  requires  certain  discontinuity  in  the  second  derivatives 
of  £  at  the  edge  of  the  tank  bottom.  For.  we  have  the  boundary  con¬ 


ditions 
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•5r~0 


V  R=l,  -U2sO 


(a) 


and 


aw 

a?  ~  bz. 


■fof  04^6.1  ,  2=-L 


(b) 


Hence 


wo^r£|,z=-l  (o 


Now  if  were  continuous  everywhere  we  must  have,  from  (a), 

that  at  R  =  („E  =  -L  •  Hence  (c)  yields  ^(^)=0 

lt  R-l,  B-~L  Ua*  "COttJt.  at  R=|;  2=-L  '  T™8 

is  Inconsistent  with  the  condition  VsJ  —  O  at  fbe  same  point.  Hence 

a1  <5 

we  mast  conclude  that  must  not  bo  continuous  at  the  corner 

£  2'"L>  '  re^uire*  that  the  series  (35)  and  its  first  and 

second  derivatives  be  uniformly  convergent  everywhere,  except  for 
,  which  loses  its  analyticity  at  the  corner  R  —  l  ,  2  =  "L* 
This  demand  singles  out  the  question  of  convergence  of  the  series  (35) 
to  the  foreground.  We  note,  in  particular,  that  any  truncated  series 
of  (35),  (taking  the  first  terms),  will  not  have  the  desired  property. 


*  The  linearised  equation  can  admit  only  one  boundary  condition.  The 
condition  is  chosen  in  Eq.  (  24-  )  on  physical  basis. 
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The  difficulty  considered  here  seems  to  be  caused  by  the 
linearisation  of  the  governing  equations.  If  linearization  was  not 
made,  the  kinematic  conditions  on  the  free  surface  and  on  the 
bottom  are,  instead  of  Eqs.  (  2 2  )  and  (  23  )•  the  following: 

J.  aH  ,  I  3H  3$  _ 

,  i_3VJ 

^•c  R'sb  tie  ~Tz. 

In  the  nonlinear  form  the  problem  is  undoubtedly  much  more  difficult, 
but  the  particular  problem  referred  to  above  disappears. 

If  we  relax  the  demand  for  an  exact  solution,  we  could  turn 
to  variational  methods.  In  the  following  the  variational  principle 
and  approximate  solutions  will  be  discussed. 


*5 
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6*  The  Variational  Principles 

If  we  multiply  Eq.  (18)  by  £<£  ,  (19)  by^TH  ,  (20)  by 
,  and  integrate  over  the  entire  fluid  volume  V  ,  the  free  surface 
S,  .  the  bottom  surface  ,  respectively,  and  transform  with 
appropriate  integration  by  parts, we  obtain 


■t^o  V,w] dss  -  ^  Iff -|*)«  **  *  £  jf 
- 1  ^  ^to)ds,  -  o 


(35) 

where  S2  is  the  rigid  side -wall  surface,  P,  is  the  edge  of  St  • 
is  the  edge  of  •  This  suggests  the  following  functional 

where 


(36a) 
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-  -+(*-+^W}W  (36c) 

^  =  4  t,  ;  *'  d2'  +  0*  4  ?Wi  6'X' (36d) 

It  is  easy  to  verify  that  the  variational  equation 

^  ^  “  °  <3?) 
undor  arbitrary  variations  of  <£>.,  Hj  W,  with  tte  stipulation  that 
vanish  at  X-Xx  ,  and  X-T-^  .  yields 

£  ! ‘ l  V’§ Si "  -  i;  l M> dS, 

M  w-WMjjswds, 

■^3 


)5Wd^-0 


(38) 
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Since  are  arbitrary  over  V,  S,,Sj  respectively, 

we  obtain  the  differential  aquations  (18),  (19),  (20),  and  the  following 
boundary  conditions: 


On 

S,  ;  Either 

or 

On 

S,  1  Either 

st  -  0 

On 

S3:  Either 

or 

-  o 

(39) 

On 

;  Either 

or 

On 

H  ;  Either 

aw  *F, 

—  -v  —L/  -  Aor 

S  w- 0  . 

The  conditions  on  the  left  hand  column  are  the  natural  boundary  condi¬ 
tions;  those  on  the  right  hand  side  are  the  rigid  boundary  conditions. 
From  physical  considerations  of  our  problem  we  impose  the  natural 
boundary  conditions  over  -S, ^  S-j  and  fj  ,  but  the  rigid  boundary 
condition  \M  =  O  over  . 

The  terms  in  the  functional  X  have  the  following  physical  sig¬ 
nificance.  The  first  term  represents  the  kinetic  energy  o:  the  fluid 
(in  dimensionless  form).  The  terms 


represent  the  change  of  surface  energy  of  the  free  surf-.te  change 
of  elastic  energy  of  the  bottom,  and  the  kinetic  energy  of  the  bottom, 
respectively. 
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To  interpret  the  moaning  of  the  other  terms,  we  write  the 
pressure  in  dimensionless  form 

P(  zX)  =  "  i  It  *  %  ^  2  +  ^  ^  J  (40) 

Then  the  work  done  by  "P  through  a  displacement  Y4  is 

^P(R,eL2,r)d2 

On  evaluating  this  integral,  and  neglecting  third  order  infinitesimals, 
we  see  that  it  is  equal  to 

[Noto  that  ^  ]>  Similarly,  the  terms 


-  4(-rC,@j..LW  '  W«£] 

are  equal  to  tho  integral  of  P(2,P.S,T)  between  ~l.-t bj  and  -  lj  , 

with  the  sign  reversed  because  the  pressure  P  acts  on  the  upper  side 

of  the  bottom.  Finally,  A— reorosents  the  work  done  by  the 
Jsi 

inertia  force  due  to  the  gravitational  acceleration  (q(T  )  through  a 
displacement  W  • 

^  are  tho  potential  of  the  vertical  forces  (positive 

upward)  acting  on  the  free  surface  and  the  tank  bottom,  at  the  edges 
P,  and  FI,  respectively.  F",  arises  from  the  capillary  surface 
energy;  pj  normally  arises  from  the  reaction  of  the  wail  on  tho 
membrane,  but  it  can  be  imposed  by  an  external  agency.  If  P,  is 
assumed  to  be  proportional  to  H  .  then  the  so-called  capillarv- 
hysteresis  is  obtained.  Eq.  (14)  and  Eq.  (25)  presuppose  F,  —  0  • 
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7.  An  Alternate  Form  of  the  Variational  principle 

It  is  Known  that  an  approoriate  functional  from  which  the  equa¬ 
tions  of  fluid  mechanics  can  be  derived  from  the  calculus  of  variations 
is  the  "pressure  energy" 

L,=-  jYdv  (4i) 

integrated  over  the  entire  volume  of  the  fluid.  See  Bateman,  (Ref.  IS) 
Wang  (Ref.|4  )•  With  this  information,  we  consider  a  functional 


I  ~  Li  -t  Li-V  L-J  (42) 

where  L,  is  given  by  (41),  with  2  expressed  in  Eq.  (40);  whereas 
and  are  the  Lagrangian  functions  of  tlie  free  surface  and  the  bottom 
membrane,  respectively: 


hi—  ftj  S  (43) 


It  can  be  verified  at  onco  that  X  so  obtained  from  (42)  is  exactly  what 
is  given  in  (36),  except  for  the  term 


which  mahemo  contribution  to  the  variational  equation  (37).  It  is  im¬ 
portant  to  note  that  the  volume  integral  in  (41)  means 


^  de 


-  tAAl 


The  variable  limits  for  0  must  bo  accounted  fori 

In  practical  applications,  it  is  common  to  choose  <&  which 
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in  Eqs.  (36),  (36a)  as 


Then  X  contains  surface  integrals  only. 
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3.  An  Approximate  Soluticn 

An  approximate  colution  based  on  the  Rayleigh-Ritz-Galerkin 
procedure  will  be  given  below.  We  chooee  §  H  and  \a/  in  the  fol¬ 
lowing  form 


$  =  d  jS)  5  i-  C<&)  1  [c*.« 


sUU- 


+  d-tT^ 


t\Lj 


J.M-) 

cHOTl 


(45a) 

(45b) 


W  -  cle(t)  -  [  cj&  -  d„trt]  X  k  e)  -v  (X + 0  ['  - 


(45c) 


which  satisfy  Eqa.  (18),  (21),  (22),  (23)  and  (25).  Eq.  (24)  is  satisfied 
by  taking 


d.co  ■=■  -  cfnCrV)  T.(k) 


(46) 


Eq.  (19)  is  satisfied  if  we  impose 

si?  cicc)  -v  B>,  6fr)  dfr)  =  6  (47a) 

-S2_^  c^r)  +  •tv.'ClVvX  ckfr)  -  0  (47b) 

Eq.  (20)  is  not  satisfied,  whereas  Sq.  (38)  now  becomes 

f  {  V*W  O 

55  (48) 

Now  we  shall  satisfy  (20)  in  the  sense  of  (48)  by  choosing 

^VM  -  CXCU -j.Ct.RjJ  <49> 


Then  (48)  yields 


ni,  ( _ ? -SLfcwi-  +xr«)  C„  -  j£  (tUlv.ux^Vdn  - 

/.  ,  ^  5lv\  d .  s-2.cC  T i(5b*')K'\  \  .  Bm_  ~1  ■) 

=•  -  1  ■*  ^H4rH(3o, 

A  substitution  of  (46)  into  (50)  gives 


ft*  (m- cw« - ^ dW) •+  k (t*  - [c^i -4<o]  (51) 

= -  /tt^_wA 


where 


A-  O&Kl*  iXA  +  fc-L- 
V*  —  fca^iKL  -^X*,,-*  KL. 


(52) 


The  central  problem  lies  in  the  solution  of  Eqs.  (47b)  and  (51). 
These  equations  can  be  put  in  a  neater  form  by  letting 


(53) 


,  C„  -  $c, 

Then  wa  have 

(54) 

where 

>M*) .  M,=  f^ 

V$)£  +^?  y 


(S4-Q,) 


H,  =  I  >C,  +•£•  \ 


)f*'\ 

Jl  J  Wuy 
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If  is  positive  definite  and  has  two  distinct  positive  eigen¬ 

values  Si}  ,  72-\  then  there  exists  a  nonsingular  matrix,  U  ,  such  that 

Si}  O 

'  -  ‘  '  (55) 


u  Ku=  U'  n,')=D- 

Let  Y-t-TXthen  Eq.  (56)  assumes  the  normal  form 


Y  -v  DY-+  <*aPCtsVuu  Y  U*'M^(^),(56) 


which  will  be  studied  in  the  following  two  sections. 
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9.  Free  Osculations 

Consider  the  special  case  of  free  oscillation  of  the  system,  i.o. 
the  case  C(-Oi  or  =  t  •  In  this  case  Eqs.  (47b),  (51)  become  two 
linear  differential  equations  with  constant  coefficients.  By  putting 

«'-Ceiat  ,  cMc^D^e^ 

we  get 

-  Ry  a*C,  •*  D«  =•  o 

{-  /xn  ni  Rl  -  Mtf-fcuTi)  c„  Et«o. 

The  eigenvalues  of  J7.  are 


It  is  recognized  that  is  the  circular  frequency  of  sloshing  of  a 

liquid  in  a  rigid  container.  In  case  ,  (which  is  usually 

the  case  for  a  rocket  in  Earth's  gravitational  field),  we  can  write 
in  the  following  form 

Bt-f-  1'-  -  1  ,S7*’ 

SutCO  iit;,)  ,  \J2,  |  is  always  less  than  ^  ;  which  means  that  the 

tank  flexibility  lowers  the  natural  sloshing  frequency.  Note  also  that, 
when  either  — »c*>,  (rigid  tank)  o’-  >  ©o  (  the  deep 

water  case,  or  a  high  mode)  tends  to  0  . 
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10.  Stability  of  the  Solution 

To  study  the  stability  of  the  solutions  of  Eq.  (56),  (Refs.  11,  12) 

$ 

we  shall  consider  tne  following  more  general  system  of  equations 


ut 


3.  *  £ 


3.  A* 


lint 


a*-° 


unt 


g, -t  jl  %,e“"VaS 

^  "»-«  1  O*  ut-iw 


(58) 


where  E\,^  As  ,  A«>  0  ,  ,  A-i  =•  &i  0  ,£P. 

(whore  =  <rr  (B„|  )  is  an  absolutely  convergent  series. 

Eqs.  (58)  are  invariant  when  T  is  changed  to  t-i  IT  ;  therefore  if  J£t) 
is  a  solution  of  Eqs.  (58),  is  also  a  solution.  By  Floquet 

theorem  (Ref.  9),  Eqs.  (58)  have  solutions  of  the  following  form 

e^  yr) 

is-  (i)  ,  t  -  (l) 

Htt)  is  periodic  function  Vttod(Tl) .  If  a"?eo  os  an  unbounded 

solution  exists,  which  is  said  to  be  unstable.  For  a  periodic  solution 
motion  to  exist,  must  be  equal  to  an  integer,  whereas  Bi  —  O. 

Let  uo  assume  a  solution  of  the  following  form 


'So 


where  21  Yf*  "*1  O"  1^,  are  absolutely  convergent  series. 
\l»-co  '  H--to 


*  This  is  an  extension  of  Hill's  method  (see  Ref.  12,  p.  413)  to  a 
system  of  two  equations. 
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Substituting  into  equations  (58),  v;e  get 


Z  +s  2  T^e  4 

VW-»c  t\»-oo 


2.  An»_,€“  2  ^ 

u*i-»o  V^r-c>c 


(2HA  * 


—  D 


2.  r^,  C^)  e  ^ 


2,  ^*’‘<'2  y>  (59) 


-v  T  3  e2"'*  T  r  efsm+^jr 


On  rearranging  the  torms  of  the  absolute  convergent  series,  and 
equating  the  coefficients  of  to  zero,  we  obtain 


w.  (ttW' 

*'*  *iM  A  "  s 


Oo 


s  -*  2 


T~1«-VV\  0 


-r. 


!«•»> 


/!.-  4** 

(*5-in)  ^  r 

"&.-<n»r  e..-^ 


(60) 


-o 


provided  that  f\ 4ai  a, o  The  divisors  Ar^j 

are  introduced  in  order  to  make  an  infinite  determinant,  which  will 

be  formed  below,  to  be  convergent. 

Equations  in  (60)  are  a  set  of  homogeneous  equations.  For  7^_ 
to  have  nontrivial  solutions  the  determinant  formed  bv  the  coefficients 
of  the  equations  must  vanish.  Call  thi3  determinant  £(.?  ) ;  then 
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\oi^ |  =-  o 


(61) 


where 


—  A.  -  (<  l 
fto  -  4-  W1 


*2 


**■  7T  "  ^  -v 
60- 


^2hvn  ~ 
C£j_fc*\,  \\  SS. 


gt- 

S»  - 

Aim-  v- 


Ao“  ■A-V"1' 

-*x 

&«  -  4  Mv 


".''-  “41,12 -  i 

We  consider  another  infinite  determinant  Zh,(.i  -  1(6; 


~(or  1<* -W  +  0 
•f <n  2  vn  ti  -  h  *■  o 

where 


£».«  -  ' 
^•v"  — 


(62a) 


Aa*.-* 


f\0  -  , 


JK.-K+-0 


(62b) 


R.  _  &mt|-n 

ri«*i ,  h  —  - —  — - — 1 — 

Px***^*^  &c(iVJMj[(  (62c) 


Since  TT  2  I  66„\  converges  provided  ^  does  not 

have  such  a  value  that  one  of  the  denominators  of  (i  ^  }  vanishes. 
Thus  the  Infinite  determinant  is  absolutely  convergent. 

Then  (Ref.  7) 

,  .  /).  ~  fA,-Cl-a«rlT 

J,1  ,»..J 

v-iO^~A'lv-I  (>T*3a.)  4-  -^6.) a— ^ 

*~W>  £~"(SK) 


(63) 
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We  note  some  interesting  properties  of  :  (1) 

is  a  meromorphic  function  of  and  tends  to  I  as  ?  i’OO; 

(2)  AjSl)  is  a  periodic  function  of  ^  with  period  zi  .  If  we  form 
another  function 

K,[^  m+fc)  -  uK  l  (cf-J As)] 

-life*  *c»'5  at  |61  -!&«,)]  <64> 

—  Effect  ^ 

where  S  are  so  chosen  that  )  ha3  no  poles  at  V?=  . 

Since  C)  is  a  periodic  function  of  ^  ,  it  follows  that  F(\)  has 
no  poles  at 

i'?  =2n  ±Ja.  ^  zn±IB„  ,  n  =  ±i,±z . 

Thua  is  a  meromorphic  function  with  no  pole  on  the  entire 

plane.  Ft'S'  is  certainly  bounded,  therefore,  by  Liouville's 
theorem.  F(F)  must  be  a  constant,  say  C  .  As 
Therefore 

c*l  +  2^x+E*V  *» 

=  1  -2  CB.tE.oi  «> 

Hence  2t*S,*0  ,  and  |  for  all  £  .  Using  this  result 

and  Eqs.  (63)  and  (64)  we  get 

A(i.^)=  -zS,1—  *Si  r^5^~ 

(?-K>)  ( ll^o) 


(65) 
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where  5j  are  some  constants  relating  to  S'/i  r  A0  and  • 
Put  i'jio,  and  1  in  (65)  and  we  get 

S,*  am 

=  \-+  MJJ*)  i  ^  - *">]  (66) 

^  ^3-  i'1'  ^ 

In  a  special  case,  if  the  coefficients  of  (58)  are  even  functions  of  t  » 
then  (58)  is  unchanged  when  we  change  't  to  - 1  ;  we  see  that  if  T 
is  a  solution,  then  -~£  is  also  a  solution.  Therefore,  if  —  o 

when  0  ,  then  ^-*.0  .  Therefore,  when  we  want  to  find 

the  roots  for  A(rt)-°  ,  we  always  have  jU*(fl.'V)so  and  the  roots  of 
(65)  can  be  written  out  in  a  simple  form. 

Is:- S’,  {67) 

For  a  bounded  solution,  i.  e.  for  (2»  ,  we  must  have 

IZ  ±  nni.  2-0  (68) 

In  Eq.  (65)  by  putting  Al»\)  — 0  ,  we  can  compute  ^  ,  and  de¬ 
termine  whether  this  is  an  unbounded  solution  or  not.  Then  from  (60) 
we  can  compute  ^  ,  and  obtain  the  complete  solution  of  Eqs.  (58). 

For  a  periodic  solution  we  must  have  ,£,16>  =  0  or  AtO^c. 

If  a  periodic  solution  of  an  inhomogeneous  counterpart  of  Eq. 

(58)  is  considered,  e.g. 

%-r  &n 


where  Btty  is  a  column  matrix  with  it3  elements  as  periodic  function. 
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we  can  use  (59)  by  putting  o  ,  whereas  Eqs.  (60)  become  inhomo¬ 
geneous.  If  ^6)-^.  o  ,  we  can  solve  for  Y^.  uniquely.  If  A(o)— e> 
we  are  on  the  boundary  where  Eqs.  (58)  have  an  unbounded  solution. 
Therefore  for  such  an  inhomogeneous  equation,  as  in  our  problem  in 
Eq.  (57),  the  zone  of  instability  is  determined  by  the  homogeneous 
solution. 
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11.  Concluding  Remarks 

In  the  above  the  formulation  of  the  problem  of  stability  of  a 
liquid  in  a  cylindrical  container  is  given.  The  analytical  difficulties 
that  occur  in  this  problem  are  pointed  out.  Variational  approach  is 
favored  for  approximate  solutions.  The  resulting  mathematical 
analysis  is  carried  out,  but  numerical  results  are  not  included  in 
this  report.  The  trends  of  fluid  stability  in  low  gravity  conditions, 
as  influenced  by  tank  flexibility  and  surface  tension,  will  be  pre¬ 
sented  in  a  future  report. 
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Fig.  1. 


Geometry  of  the  Problem  in  Nondm.ensional  Variable. 


